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Abstract: The drag of a heavy quark and the jet quenching parameter are studied in 
the strongly coupled J\f = 2* plasma using the AdS/CFT correspondence. Both increase 
in units of the spatial string tension as the theory departs from conformal invariance. 
The description of heavy quark dynamics using a Langevin equation is also considered. 
It is found that the difference between the velocity dependent factors of the transverse 
and longitudinal momentum broadening of the quark admits an interpretation in terms of 
relativistic effects, so the distribution is spherical in the quark rest frame. When conformal 
invariance is broken there is a broadening of the longitudinal momentum distribution. This 
effect may be useful in understanding the jet distribution observed in experiments. 
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1. Motivation 

Experiments of heavy ion collisions at RHIC show a quantitative deviation from pertur- 
bative QCD predictions of heavy quark jet production [1]. There are several indications 
that a deconfined state of matter that behaves as a strongly coupled plasma is formed at 
the center of the collisions. For instance, hydrodynamic simulations using a very low vis- 
cosity have been quite successful in describing the observed elliptic flow at low transverse 
momentum [2]. Another example is the disappearance of back-to-back jets, that could 
be understood as the effect of a highly dissipative medium. Further evidence comes from 
lattice theory [3], where it was observed that at the energy densities reached at RHIC the 
equation of state of QCD is still far from the free gas value. 

Some of these observations agree with predictions of AdS/CFT duality [4] that pro- 
vides a holographic description of some strongly coupled gauge theories. For instance, the 
pressure of the strongly coupled N = 4 super Yang-Mills theory was predicted to be about 
75% of the free gas [5], which is similar to the results found in the lattice theory in the 
range of energies of RHIC. Another prediction was that the plasma shows a hydrodynamic 
behavior with a very small viscosity [6]. In the context of heavy quarks, a classical string 
computation in AdS / CFT gives predictions for the drag force [7-9] , and there are different 
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estimates for the energy loss of a fast particle moving through the medium. The various 
approaches involve studying string configurations corresponding to massless quarks [10,11], 
the jet quenching parameter from a dipole source [12] or the correlation functions of charged 
currents in the plasma [13]. The behavior of the drag and jet quenching in a non-conformal 
theory will be the main interest of this paper. 

Both the drag force and the jet quenching have been studied in a large number of 
situations [14-40], including some non-conformal theories. As opposed to the shear vis- 
cosity, these quantities do not show an universal behavior. Although according to the 
lattice results the equation of state of QCD in the RHIC range is not very far from a 
conformal theory, it is necessary to quantify how the deviation from conformal invariance 
will affect the energy loss of the jets. For this, one needs to use a framework that can 
connect smoothly with the conformal case. There have been several recent attempts in 
the context of AdS/QCD that show an interesting behavior [38-40] , but it is certainly 
desirable to study a case where the deformation is fully understood. For this purpose, the 
J\f = 2* theory [41,42] is an ideal laboratory. Conformal invariance is explicitly broken by 
the introduction of a mass deformation in the N = 4 theory. The ultraviolet fixed point is 
the N = 4 theory, so the holographic correspondence is very well understood. The ther- 
modynamic analysis [46, 47] shows that even when the equation of state starts deviating 
significantly from the conformal case there are no phase transitions. Therefore, one can 
tune the mass from zero to quite a large value and study the variation of the observables 
of interest. 

Although the numerical results are given for the N = 2* theory, the analytic results 
presented here do not depend on its particular properties, and should be general for a large 
class of models that describe scalar relevant deformations of a strongly coupled conformal 
theory in four dimensions. One should add that the deformation must have a dual descrip- 
tion in terms of a scalar supergravity field, which it is not necessarily true for all of the 
possible relevant deformations in the field theory. 

The paper is organized as follows. In section 2.1 there is a small review of the main 



features of the N = 2* holographic dual. In section |2.2| some general properties of the black 
hole solutions and the thermodynamics of the N = 2* theory are discussed. The choice of 



string configurations used to compute the different quantities is explained in section ^3 



Section || contains a derivation of the spatial string tension that will be used as reference 
scale. The next sections contain the computation of the drag coefficient and momentum 
broadening (section^), and the jet quenching parameter (section |5|) . The numerical method 
used to calculate the values in the Af = 2* case is explained in the appendix |A[ Results 
are summarized and discussed in section H. 



2. Preliminaries 

2.1 The supergravity dual of N = 2* 

The contents of this section can be found in the references [42-44], but are included here 
for completeness. The field content of N = 4 can be grouped in J\f = 1 superfields as a 
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vector multiplet and three chiral superfields with scalar components 



&i = X 2i - 1 +iX2i, i = 1,2,3. (2.1) 

In the J\f = 2* theory, two of the chiral fields 3>i and $2 a re grouped in a N = 2 hypermul- 
tiplet and a mass is introduced in the Lagrangian. The remaining fields form an J\f = 2 
vector multiplet. The mass is introduced through the following relevant operators: 

2 = tr (-X 2 - X\ - XI - X\ + 2X1 + 2X1), 

03 = tr (A1A1 + A2A2) + (scalar trilinear) + h.c. . 

where Aj are the fermionic components of the multiplets. This choice of operators corre- 
spond to the SO(6) representations 20' and 10 + 10, that map to two different scalar fields 
in the holographic dual, a and %. The last can couple to the scalar singlet mass operator, 

Q 1 = tr (Xj + X 2 2 + X 2 + Xl + Xl + Xl) , (2.3) 

that is not BPS protected and has no associated supergravity field. The mass term comes 
from the combination 

Oi ~\°* = \ tr ( X i + X l + *3 + Xl) • (2.4) 

The dual theory is described by solutions of five-dimensional N = 8 supergravity with 
the scalar fields turned on. This is a consistent truncation of type IIB supergravity in ten 
dimensions, and it is possible to uplift the five-dimensional solutions to ten dimensions, 
although this is far from trivial. In the zero temperature case the five-dimensional metric 
can be put in the form 

dsl A = e 2A{r) {-dt 2 + dx 2 ) + dr 2 . (2.5) 

Thanks to supersymmetry, it is possible to obtain the function A(r) and the scalar fields 
a(r), x(r) as solutions of a system of first order differential equations. 

dr 3 ' dr 4 da ' dr 4 d X ' 



where 



W = - e - 2a/V * — cosh(2x) • (2.7) 

is the superpotential and g is the five-dimensional gauge coupling. The value of g sets 
the radius of curvature L = 2/g, related to the 't Hooft coupling in the dual theory via 
L 4 = g 2 YM N(a') 2 = X(a') 2 . 

The zero temperature theory has an interesting structure that in the gauge theory can 
be formulated in terms of the Coulomb branch of the moduli space. On the gravity side 
the supergravity solution has a singularity at the origin of the space. This singularity is of 
the good kind [45] and from an analysis of D-brane and string probes in the geometry, it 
is possible to relate it to a particular configuration in the Coulomb branch of the J\f = 2* 
theory where all dyonic states become massless on a ring (the enhangon) around the origin 
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of moduli space [43,44]. At finite temperature the Coulomb branch is lifted and there is 
no enhangon singularity in the geometry. 

At non-zero temperature supersymmetry is broken and it is necessary to use the second 
order equations of motion [46,47]. 

ldV ^ IdV 



The potential for the scalar fields is 



V = 9 - 
16 



dW\ 2 fdW x ~ 

+ 



da J \ dx 



a 2 

^-W 2 . (2.9) 



The metric can be written as 

dsl, = e 2A ^ {-e 2B ^dt 2 + dx 2 ) + dr 2 , (2.10) 
and the functions A(r) and B(r) can be obtained from the Einstein equations 

= - \ 9iXV TV (2.11) 

where 

V = 3^ad v a + d^xduX ~ [d^a + d pX d p x + V] . (2.12) 

is the energy-momentum tensor of the scalar fields. 

From the equations of motion one can deduce the asymptotic behavior of the scalar 
fields. As r — ► oo 

x = ke-rl\l + ...), p = = 1 - \k 2 \^ +■■-, (2.13) 

where k = mL is the mass of the hypermultiplets. This radial dependence corresponds to 
scalar fields with masses 

m 2 x L 2 = -3, m 2 a L 2 = -4. (2.14) 

This matches the relation between the mass of a scalar field and the conformal dimension 
of the dual operator, m 2 L 2 = A(A — 4), with A = 3 for the O3 operator and A = 2 for 
the O2 operator. At the horizon r — > the scalar fields have a constant value 

X = Xo + Xir 2 + • • • , p = Po + Pir 2 + ... . (2.15) 

As was mentioned before, it is possible to uplift the five-dimensional solutions to a full 
ten-dimensional geometry. In the Einstein frame it looks 

ds 2 10 = n 2 dsl yl + ds 2 5 , ( 2 - 16 ) 

where 

ds 2 = ^ (c-V + f cos 2 e(^L + ^Ipl) + sin^ 8*£) (2.17) 

p \ \CX2 Ai J X 2 J 
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with the warp factor 0, c and p defined as 

~ 2 = {cX x X 2 )V^ c = cosh(2x)i p = e */V3. (218) 
P 

The scalar functions are 



Xi = cos 2 9 + p 6 csin 2 
X2 = c cos 2 9 + /) 6 sin 



2.^6^2. (2-19) 




The one-forms cij are the SU(2) left-invariant forms satisfying dai = e ijkCj A cifc. 

For the computation of string solutions, one also has to take into account the dilaton 

(2.20) 

The axion and the three and five form fields are also present in the geometry, but they will 
not be important for the analysis presented here. This is the case because the components 
of the NS flux along the directions of the string probes studied here vanish. In other cases 
the situation can be different, an example of this would be the non-relativistic backgrounds 
of [48] . This would affect for instance the considerations of section [2.3j . 

2.2 Black hole properties and thermodynamics 

Some general properties will be presented here. For concrete calculations the equations are 
solved numerically. Details are explained in the appendix [A]. 

First recall the form of the five-dimensional black hole metric 

dsl ! = e 2A ^(-e 2B ^dt 2 + dx 2 ) + dr 2 . (2.21) 

For asymptotically AdS spaces, the behavior as r — > 00 is 

A( r )->j+Ao, (2.22) 

and 

5(r)->0. (2.23) 

The constant Aq is set to zero by a rescaling of the spacetime coordinates. It is always 
possible to choose a coordinate system where locally the metric component g rr = 1, so 
this form of the metric is completely general for a geometry that has spatial rotational 
invar iance. 

The horizon is at r = 0, where the blackening function e 2B vanishes. For a regular 
horizon the metric functions have the following expansions 

eB(r) = 6o z( 1 + 6 4 + -) ' (2 ' 24) 



e^W =a0 (i + Ol l_ + ... j . (2.25) 



2 
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The smoothness of the Euclidean solution and the area of the horizon determine the tem- 
perature and entropy density 

T = m> s = ^- < 226 > 

(G5 is the 5d Newton constant). Although in the uplift to the ten dimensional metric ( 2.1€| ) 
the internal space and the conformal factor depend on the scalars, one can show that they 
do not affect the values of T or 5 1 . 

In AdS§ {N = 4 theory), the values are 

ao = ^, 60 = 2, (2.27) 

where uh is the horizon radius in Schwarzschild coordinates and 6q = 4 is related to 
exponent in the blackening function /(it) = 1 — ujj/u 4 , that in the AdS space is fixed by 
the number of dimensions. 

In the general case, the entropy can be cast in a more useful form in terms of the 
N = 4 entropy. Using 

L 3 N 2 , , 

ift < 2 - 28 > 

the entropy density is 

S = ^N^ = y M= ,. (2.29) 

So the ratio 2/6o controls the deviation from the conformal case. 

One of the Einstein equations can be integrated, giving a relation that will be useful 
for the computation of the drag coefficient. Notice that for scalar fields that depend only 
on the radial coordinate, the combination 

QuRoq - OooRa = 0, i = 1, 2 or 3 , (2.30) 

is independent of the scalar fields. This gives the equation 

B" + (4A' + B')B' = => In B' + AA + B = const . (2.31) 

This formula is true for any Lorentz-invariant deformation involving scalar fields, not only 
for the N =2* theory. 

These formulas are quite general, but in order to describe the thermodynamic behavior 
of the J\f = 2* plasma one has to solve the equations of motion and find the right geometries. 
This was done in [46,47] and later on, hydrodynamic properties of the N = 2* theory were 
also studied in [49-52]. The free energy was computed numerically up to values m/T ~ 7. 
A good approximation is given by the fit 

Fm=2* ^ F M=A e~ m /W (2.32) 

For a value m/T ~ 6, the speed of sound was found to be vl ~ 0.825 fj^pj., with Vq FT = 
1/3. This implies that the equation of state of the N = 2* plasma is quite close to the 



1 The area of the horizon just gets a factor that is the volume of the undeformed S 5 . 
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conformal case up to values of the mass m/T ~ 7. If the bulk viscosity over shear viscosity 
ratio C/rj approximately saturates the bound observed in [51] 



c^A ..2 



then, at m/T ~ 6, Q/r] ~ 0.12 is quite small, but starts to be comparable to the shear 
viscosity. In the numerical evaluation of quantities, the mass will be taken up to larger 
values m/T ~ 17, so the theory is in the region where deviations from conformality should 
be substantial. 

2.3 Averaged Wilson loops 

The analysis of the energy loss of heavy quarks in the plasma involves the computation 
of some Wilson loop configurations. In = 4 it is possible to use a classical string 
configuration in the holographic dual to evaluate the expectation value of a generalized 
version of a Wilson loop, involving the gauge and the scalar fields [53] 

(W(C,6)) = Ulc d ^ A ^+ eIx ^)\ . (2.34) 



where x^(r) parameterizes the contour C of the Wilson loop and the 9 1 I = 1, . . . , 6 are 
components of a unit vector in R 6 , so they parameterize a five-sphere S 5 . In the holographic 
dual the S 5 is realized geometrically, and the generalized Wilson loops can be described 
using classical string configurations that sit at a point of the S 5 . In the TV = 4 theory the 
SO(6) R-symmetry that acts on the S 5 is unbroken, so the values of the Wilson loops are 
independent of the value of 9 1 . 

In the M = 2* theory the moduli space is partially lifted by the masses, so the S 5 in the 
holographic dual is deformed. The SO(6) ~ SU(4) symmetry is broken to a SU(2) x U(l) 
subgroup. In this case one should be careful with the choice of Wilson loop, since the 
expectation values now depend on the values of 9 1 . Other situations where this can happen 
include the R-charged black holes, as shown in [54]. A possible definition of a Wilson loop 
that is independent of the 9 J 's is an averaged Wilson loop 



(W(C)) = dn 5 (9) (W(C, 9)) . (2.35) 

Where the average is over the parameters 6 1 that define the Wilson loop in the field theory. 
This corresponds to an object with zero R-charge, so it is closer to the properties of a Wilson 
loop involving only gauge fields. A similar kind of average was considered in [55] for Wilson 
lines in the N = 4 theory with dependence on the internal directions, as a way to estimate 
the quark-antiquark potential for QCD. 

Clearly in the N = 4 theory the averaged Wilson loop has the same expectation value 



as (2.34). Although it is possible to find an explicit expression for ( W(C)), one can use a 



saddle point approximation instead. For this, one writes the Wilson loop as 

(W(C,6)) = e iS{C ' e) or (W(C, 9)) = e- s{efi) , (2.36) 
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depending on whether the Wilson loop is timelike or spacelike. The action S(C, 8) computed 
in holographic models has a large y/\ factor, so in the strong coupling limit the integral 
over the S 5 volume will be dominated by minimal action configurations 2 

(W(C)) ~ (W(C,6 mhl )) . (2.37) 

where 6 m [ n is the value that minimizes the action. There is actually a full set of values 
related by the remaining SU(2) x U(l) symmetry, but as for the N = 4 case, it is enough 
to take one representative. Notice that this formula would be valid only as long as the 
string configuration is at a fixed point in the compact space or moves along an isometric 
direction. That would not be the case if for instance an NS flux induces a force on the 



string, as was mentioned in section 2.1. The minimal action configurations at a fixed point 
will be the ones used to compute the string tension, drag and jet quenching. 

A generalized Wilson loop ( p. 34 ) has a holographic description in the semiclassical limit 



as a classical string configuration in ten dimensions. The string is a two-dimensional surface 
with a boundary at r — > oo that is the contour of the Wilson loop C. The configuration is 
determined by the equations derived from the Nambu-Goto action in the Einstein frame, 

SNG= 2~hJ d2(jC = 2~h/ J d2<J ^ ~ 6lP/2 ^ 2 ' (2 ' 38) 



where gi is the pullback of the five-dimensional metric ( |2. 10 ), ip is the dilaton ( |2.20| ) and 



ft is the warp factor From (|2.l3 ), the factor Q 2 -> 1 as r — > oo. On the other hand, 



when r — > the scalar fields tend to a constant ( [2,15 ) and the factor Q 2 in the string action 



depends on the position of the string in the internal space. 

String configurations with minimal action are localized at 9 = ±7r/2, that corresponds 
to the part of the geometry associated to the M = 2 moduli space. In this case the factor 
is minimized for all values of r, 



2cosh(2 x) 
V l + cosh 2 (2 X ) 

So there will be a dependence on the value of the scalar fields. In the r — > limit, the 
factor becomes a constant: 

2 _ 2 / 2cosh 2 (2 X0 ) 
^- P °Vl + cosh 2 (2 Xo) - (2 - 40) 



3. Spatial string tension 

Spatial Wilson loops in the deconfined phase show an area law behavior when the size of the 
Euclidean time is very small compared to the size of the loop. It is interesting to compute 
the spatial string tension using holography for two different reasons. In the first place, it is 

2 Strictly speaking, for timelike Wilson loops the condition is that the action is extremized, but in the 
cases studied here the physically sensible saddle point corresponds to the position on the S where the 
action is minimized. 
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a physical quantity that one may compare directly with results from lattice computations 
[3,56] (see [57] for a comparison with different methods). In a phenomenological application 
to the QCD quark-gluon plasma, one could fix the string tensions of lattice QCD and of 
the holographic model to be equal and extract the values of other parameters, in the 
spirit of previous works, as for example in [8, 19]. The second reason, and the one that 
motivates its calculation here, is that it has the same parametric dependence on the 't 
Hooft coupling as the drag force coefficient and the jet quenching, so it is canceled out in 
the ratio. In principle one could accomplish the same by taking the ratio with the values of 
the conformal theory, but notice that this does not give useful information about intrinsic 
properties of the theory. In the next sections the values of the drag coefficient and jet 
quenching parameter will be expressed in terms of the string tension. 

A spatial Wilson loop can be described holographically by a classical string configura- 
tion ending on the appropriate contour at the boundary. The contour of the Wilson loop 
is chosen to be rectangular in the x and y directions. The length in the y direction L y is 
taken to be very large, so it is a good approximation to neglect any dependence on the y 
coordinate. The choice of worldsheet coordinates is y = r and r = a, and the profile of the 
string is determined by a function x(r). The string worldsheet Lagrangian is 



C = n 2 e A Vl + e 2A x' 2 . (3.1) 
From this, the following equation of motion is obtained 

x = , . (3.2) 

Choosing the integration constant c = f2 2 (ro)e 2j4 ( r °) implies that the string profile ends at 
tq. The length in the x direction and the regularized action of the Wilson loop 3 are 

L x f°° , ce~ A 
dr 



S L.,, 



V <7T7TT -i ' 

J3L( / dr n 2 e A ■ \ - 1 - / drn 2 e A ) . (3.3) 
2 2W \J ro lV&e 4A -c 2 \ Jo J 

The length grows as tq — > 0. In that regime, the contributions to the integrals from 
r — > oo are small, since the integrands decay exponentially and the integrand has a square 
root divergence at ro- One can use the approximation e" 4 ^ ~ e A ( r °\ Defining A(r) = 
A(r)-A(r )and^=gg J , 

l, i r°° . i 



-f 

«0 Jr, 



dr 



Ly /■« 

' ir -»0 



dr = . (3.4) 



2 27TQ' J rQ ^ ^4 g 4A(r) _ J ' 



^Obtained by subtracting the straight configuration x = 0. 
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Figure 1: Spatial string tension of AT = 2* theory for different values of the mass deformation. 
The dashed line corresponds to the approximate formula ( |3.6| ) and the solid line to a numerical 
evaluation using the string solution. 



Remember that ao = e A ^\ as defined in (2.25). At large distances L x 3> 1/T the Wilson 
loop follows an area law 



s ~ a 



-L X Ly — I — 



'27TQ' " y \b J u 2 
In terms of the M = 4 theory result, the string tension is 



LxLy = CT S L X Ly 



O-.s 



ni 



(3.5) 



(3.6) 



In figure [l| this approximate formula in the M = 2* theory is compared with a direct 
evaluation using the string solutions for several values of m/T. The agreement is good 
only for small masses, although it captures the right tendency of the string tension to 
decrease as the mass increases. 



4. Drag force and Langevin equation 

In the formation of the quark-gluon plasma, soft light quarks and gluons thermalize very 
fast but heavy quarks remain out of equilibrium for a relatively much longer time if their 
mass is much larger than the temperature M/T 3> 1. At weak coupling the quarks lose 
energy through bremsstrahlung and collision processes, that can be described within the 
framework of kinetic theory. The momentum change of the heavy quark may be described 
by a Langevin equation. A quark with spatial momentum pi moving through the plasma 
experiences a force [58] 

^- = -VDPi + (Ht)+Cht). (4.1) 
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This is a phenomenological description where the effect of the medium has been separated 
in two parts, the drag ~ rjjj and a stochastic contribution ( T , ( L . The drag force opposes to 
the movement of the quark —r]DPi, where t)d is the drag coefficient. The term coming from 
random collisions with the components of the plasma & and (T broadens the momentum 
distribution in the longitudinal and transverse directions. Usually the random components 
are taken as Gaussian white noise, but in general the momentum transfer coefficients kl 
and kt can depend on the momentum 

<CP(*)C/(*0> = P -f^L(p)5(t - t>), (Cf (t)Cj(i')) = (*i - P ~f) «r(p)*(t - 0- (4-2) 

In principle, the drag coefficient can also depend on the momentum ijd = r]£>(p). The 
fluctuation-dissipation theorem relates the zero-momentum values of the drag coefficient 
and the momentum transfer kl(0) = kt(0) = k, giving the Einstein relation 

Vn(0) = ^ f . (4.3) 

The momentum transfer also determines the diffusion of the quark through the plasma. 
Consider a quark at rest at t = and x = 0. The mean squared position of the quark at 
later times will be 

2T 2 

(xi{t) Xj {t)) = 2Dt5ij, D = , (4.4) 

where D is the diffusion coefficient. 

At strong coupling there is no well-defined kinetic description of the thermal medium, 
so it is not clear that a description in terms of the Langevin equation should work. Re- 
markably, in the holographic description of a heavy quark in the M = 4 plasma, a drag 



force of the form appearing in (4.1) with a constant drag coefficient was found [7-9]. Later 
a computation of the momentum transfer showed that it depends on the momentum and 
that the Einstein relation ( |4.3| ) holds in the limit of low velocities [59,60]. 

4.1 Holographic computation of the drag 

A heavy quark moving through the plasma is described by a trailing string configuration 
in the black hole background [7-9]. The profile moves at constant velocity v in the x 
direction and it is extended from the boundary to the horizon in the radial direction r and 
from x = vt to x — > — oo in the spatial direction. Identifying the time t and r with the 
worldsheet coordinates, the profile is described by a function 

x(r,t) =vt + £(r), (4.5) 

and by the position in the internal space in the ten-dimensional geometry. As was argued 
in section a reasonable choice is to pick a point that maps into the moduli space of the 
J\f = 2 theory. Introducing ([4.5|) in the Nambu-Goto action, the Lagrangian for £(r) is 



C = -V=£ = -n 2 e A {e 2B -v 2 + e U+*Bg2f/* . (4 . 6) 

The drag force in the x direction is given by the momentum flow along the string. 
Labeling p x the momentum in the x direction, the components of the background 
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metric and g a b the components of the induced metric on the worldsheet, the momentum 
flow is 

^£ = J_^ Gxxg rr d { t) = 1 n 4 e 4A+2B ? = . (4.7) 

Where is the canonical momentum 

7Tf = = . (4.8) 

? <9£' ( e 2B _ y 2 _)_ e 2A+2B ( t/2-) 1 / 2 

From this expression one finds the following equation for the profile 



In order to have a profile that can go all the way to the horizon, it is necessary that 
numerator and denominator flip sign at the same point. This imposes a condition on the 
canonical momentum 



7T 



^ A+2B \ 2=p2B = v 2 ^\ 2=p2B . (4.10) 



Using ( 4.10 ) in (fO|), one can deduce the value of the drag force. For v > 0, tt^ > 0. 

dp x vQ 2 e 2A 



dt 2vra' 2ira' 
The expression fl4.11| ) leads to the known result in the M = 4 theory, 



(4.11) 



dp x ttVA 2 v tt/\ 2 p x _ Px 

~dT = = "^~ T M = ~ m ='M ' (412) 

where a factor of the mass is extracted from the drag coefficient /i = t]dM for convenience. 
Notice that //jv"=4 is independent of the velocity. In the general case this is not true, and 
it is interesting to extract the velocity dependence. 

In terms of the functions appearing in the metric, the drag coefficient becomes 

Mdpx = 1 2B) l/2 2 2A (413) 

At large velocities v — ► 1, the value of the drag coefficient depends on the asymptotic 
behavior of the solution r — > oo. The scalar fields do not affect to the result, since n 2 i. 
The leading terms in the metric are 



e 2B 



1 - /2?e- 4r/L + . . . , e A = e r/L (l + ...). (4.14) 



where the constant \±\ corresponds precisely to the value of the drag coefficient up to a' 
factors. One can extract the value of fi 2 using the relation ( |2.31| ). Comparing the r — > oo 
(4.14) and r — ► ( |2.25 ) limits, this gives 

ln(2/i 2 ) =ln6 + 41na , (4.15) 
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so the drag coefficient for ultrarelativistic quarks is 

*W = o^.„n2 ■ ( 4 - 16 ) 



Contrary to previous claims the drag coefficient is sensitive to the infrared physics even 
at large velocities. Technically, the high- velocity drag jl\ appears in the expansion of the 



metric (4.14p as a normalizable term so it has a similar status to an expectation value, it 



should be possible to relate it to transport coefficients of the plasma. 

At small velocities v — > 0, the drag coefficient is determined by the expansions at r — > 



(12.251) 



The ratio depends on bo and the value of the scalar fields at the horizon 



*ng. (4.18) 



v /^-kL/ bo 

The drag coefficient can be written in terms of the M = 4 coefficient: 



2(27ra') 2 ^ 

hence 



= oPo! L .mI 3 = f r-l ^=4 , (4-19) 



/4-*o = (r) ^o^=4- (4-20) 



Using (|3.6|) , the drag coefficient would be 

fj, v -+ = a s , n v -^i = \l ^rjffs. (4-21) 




For = 4, this implies that fj, = cr s , as was pointed out in [65]. However, this relation does 
not seem to hold in the non-conformal cases, because of the deviation of the string tension 
from the analytic result. The numerical value of the drag coefficient in string tension units 
as a function of the velocity for several mass deformations has been plotted in figure |2|. 
Formula ( [4.18 ) is also compared with the numerical results, with very good agreement. 



The main effect of breaking conformal invariance is that the drag coefficient depends on 
the velocity, in the N = 2* theory it is larger at higher velocities. Similar behavior was 
observed in other non-conformal theories [15]. 

4.2 Holographic computation of momentum broadening 

In the field theory side, the stochastic forces (J , C,f in the Langevin equation fl4.1| ) cor- 
respond to operators built with heavy quark fields, so it is natural to identify them with 
small fluctuations of the string due to thermal radiation coming from the black hole [59,60]. 
Notice that a force has to be applied to the quark to make it move at constant velocity, so 
the quark is not in thermal equilibrium with the plasma. However, the state described by 
the trailing string is a steady state with a constant momentum flow. 
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Figure 2: (A) Drag coefficient in string tension units as a function of the rapidity r\ = arctanh-u, 
for m/T = 0, 2.08, 3.42, 5.00, 6.94, 9.40, 12.58, 16.80 from bottom to top. (B) Ratio between the low 
and High velocity drag coefficients, as a function of m/T. The solid line corresponds to the direct 



evaluation using the string solution and the dashed line to the formula (4.18). 



In the derivation of the momentum transfer, the worldsheet black hole is proposed to be 
a holographic description of the Wilson line spanned by the heavy quark, so the steady state 
is described effectively as a thermal state of a one-dimensional theory with temperature 
T s t. The worldsheet temperature depends on the velocity of the quark and it is in general 
different from the black hole temperature. At zero velocity the quark is at equilibrium with 
the plasma, so T s t(v = 0) = T. Using a path integral approach, the degrees of freedom 
beyond the worldsheet horizon can be integrated out. The effect is to introduce a random 
force for the string endpoint at the worldsheet horizon that propagates to the boundary 
and introduces the stochastic forces of the Langevin description. This is equivalent to 
consider the Hawking radiation in the worldsheet theory. This introduces the picture of an 
ensemble of strings with fluctuating velocities. The large time average configuration is the 
classical string solution and the mean deviation is related to the magnitude of the random 
forces [61-64]. 

Using this approach, the value of the momentum transfer coefficient can be found from 
a Kubo formula involving the retarded Green's function of fluctuations of the string profile 

K a = - lim — ImG ( B fl) y, a = L,T . (4.22) 

Consider small perturbations around the trailing string 

xi = vt + £(r) + 5xi(t, r) X2 = Sx2(t,r), X3 = 5x3(t,r). (4.23) 

Expanding to quadratic order, the Lagrangian for the fluctuations is 

£( 2 ) = -g^daSxxdpdx! - Gfd a 5xidp5xi , (4.24) 

i=l,2 

4 There are also first order terms but those are total derivatives, so they do not affect to the equations 
of motion and will be ignored here. 
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where 

Qf = \fa^hh^ , (4.25) 

and 

fx = &e 2A , (4.26) 

4 e 4A+2B _ n 2 

fL = n 2 e 2A (e 2B -v 2 ) ' (4 ' 27) 

The ratio Jt/Jl f°r large values of r becomes the right one for the zero temperature theory, 
fx = (1 — v 2 )/l- The velocity dependent factor corresponds to the anisotropy between the 
transverse and longitudinal fluctuations at short times or large frequencies. However, the 
small frequency behavior should be related to values of r close to the worldsheet horizon 
ru, determined by e 2B = v 2 . Expanding close to the horizon and using ( |4.10| ), one finds 

fUn) = Mr*) (i + 2^M + , «<r„) \ 

V B'(r H ) VL(r H )B'(r H ) J 

In the N = 4 theory this gives 

h{r H )= (1 _ l 2)3/2 , Mm) = (1 _ l y2)l/2 , (4.29) 

so the zero temperature anisotropy persists at the worldsheet horizon. 
The metric h a p can be diagonalized 5 h a p = diag (—a, I /a) with 



a 



( e 2B _ v 2^ n A e 4A+2B _ n T 



1/2 



Q 2 e A + B 

Close to the horizon, the function a ~ ao(r — m) vanishes 

a ~ ao(r — rjj) + • • • = 



(4.30) 



B'{r H ) ( 2 ^M + 2A'(r H ) + B'(r H )j 



1/2 

(r - r H ) + ■ ■ ■ . (4.31) 



Requiring that the Euclidean solution is regular at the horizon implies that the temperature 
of the worldsheet black hole is T s t = olq/ (47r) . The temperature at v = coincides with the 
black hole temperature T s t = T. In the N = 4 theory, the worldsheet temperature has a 
simple dependence with the velocity T s t = T/^/7. The deviation from the conformal case 
has been represented in figure ||. 

It is convenient to use a Fourier decomposition of the fluctuations. A Fourier mode 
5x = e~ lujt cp(r) obeys the following equation (a = T, L) 

uj 2 

ip'l + (log(a/ a ))Va + — ip a = . (4.32) 
or 



5 This is done by changing the time variable t to t = t+ g{r), with g' = —h% r jh%f After this change of 
variables, the new components h a fi of the metric are h t t = htt, h tr = 0, h rr — h rr — h%,/h t t- In the text 
the new variables are relabeled to t — * t, h — > h. 
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Figure 3: Temperature of the worldsheet black hole as a function of the rapidity r\ = arctanhi;, 
for m/T = 0, 1.28, 2.08, 3.42, 5.00, 6.94, 9.40, 12.58, 16.80 from top to bottom. 



In order to compute a retarded correlator, the solution must be ingoing close to the horizon 

<Pa = (r- r H Y™ /ao F a (r) F a (r H ) = const. . (4.33) 

One can use a low frequency expansion 

F = (r) + wF a (1) (r) + • • • . (4.34) 

The leading term is just a constant and can be fixed to Fa (r) = 1. The equation for the 
next term can be simplified by choosing W a = —iF a 

W' a + (log(af a ))'W a + ' 



ao(r - r H ) \r 



One can rewrite this as 



[afaW a ]' - — 

a 



r — th 

Then, imposing regularity at the horizon r — > rjj, the solution is 



(log(af a ))') =0. (4.35) 

- r H ) 

(4.36) 



0. 



f H 1 



(4.37) 



a/a a (r - r H ) ' 
where fj* = f a {r H ). 

In order to compute the Green's function a cutoff ta 3> L, rn is introduced and the 
normalization of the solutions is fixed at the cutoff Y a {oj,rjs) = e~ tujt , so 



-icot Va{r) 



(4.38) 



^a(^A) 

The retarded Green's function is proportional to the boundary action as the cutoff is taken 
to infinity. To leading order in oj, the result is 

1 , fH 

2ira' 



GrH") = 7^1 lim fa^hh rr Y:(oj,r A )d r Y a (u J ,r A ) = -iu;^ 7 + 0(u J 2 ). (4.39) 
lira r A ^oo >«r™i 
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Using ( 4.28 ) and ( 4.39| ) in ( 4.22; ), the value of the transverse momentum transfer is 



k t = 2~fT st p . (4.40) 
The longitudinal momentum transfer can be simplified using that v = e B( - TH) and 



dp _ (dp/dr 



dv (de B /dr) 



r=r H 



il(r H )B'(r H ) B'(r H ) 



The complicated dependence on the derivatives of the functions of the metric is the same 



as in (4.28). This allows to write the value of the longitudinal momentum transfer as 



k l = 2 7 3 T st p + ^T st v^ . (4.42) 

There is an additional term due to the velocity dependence of the drag in the non-conformal 
theory, and the velocity dependent factor in the conformal theory is 7 3 instead of the factor 
7 found for the transverse momentum transfer. In the following an interpretation for this 
difference is given. 

Assume that it is possible to describe an ensemble of heavy quarks using a Liouville 
equation. The probability density at time t of finding a heavy quark in the position x and 
with transverse and longitudinal momenta pr and pl is p(x,pT,PL',t) 6 - The continuity 
equation is 

dp d , . d , d , , A 

-£ + 7T (xp) + — (p T p) ~\~ tj (plp) = . (4.43) 

at OX OPT OPL 

The trailing string describes a quark moving through a viscous medium and subject to an 
external force in the longitudinal direction F v = jvp 

PL = F v -p^, VT = -p V f l . (4.44) 

The effect of the stochastic forces can be taken into account adding to the right-hand side 
of ( 4.43| ) terms of the form ^d 2 p/dp 2 . In the absence of forces F v = and p = 0, there 
is a homogeneous solution p ~ (Kt)^ 3 ^ 2 e~ p2 ^ 2Kt ^ to (4.43) that describes correctly the 



spreading of the quarks due to Brownian motion. 

The question is now if there is a static and homogeneous solution to (4.43). This will 
mean that the heavy quarks and the plasma are in a steady state configuration with a 
constant flux of momentum, introduced in the system by the external force acting on the 
quarks. When there is no force the average velocity of the quarks vanish and the system 
is at thermal equilibrium. Under these conditions, the equations reduce to 

p T kt dp 
^M P+ Tdp^ = °> 

PL . K L dp 

-™ + i*M p + Tdj£ = ' (445) 



°AU these quantities are vectors but arrows are dropped for notational simplicity. 
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There is indeed a solution in the form of a Boltzmann distribution centered at non-zero 
longitudinal momentum 



Peq (PT, PL ) oc exp exp . (4.46) 

Focusing on the conformal case, the factor 7 2 can be understood as a consequence of the 
underlying Lorentz symmetry. By doing a Fourier transformation, the resulting distribution 
in the spatial coordinates has widths 

< A 4>~^T7^ 2 , <A*i> ~ ^—f = \l? . (4.47) 
"fTstM TT st M Y 

The square roots of the widths (Ax 2 ^ i) 1 ^ 2 give the mean length of the distribution along the 
transverse and the longitudinal directions in real space, so under a longitudinal boost the 
transverse width should not change and the (square root) of the longitudinal width should 
suffer the usual Lorentz contraction. If the distribution of the ensemble of heavy quarks 
is spherical with mean radius L in the rest frame defined by the average values (xt) = 0, 
(xl) = vt, then the boost to the rest frame of the plasma will account for the different 
7 factors in the longitudinal and transverse widths in ( 4.47| ). The possible relation of the 



additional j 2 factor in the longitudinal direction with Lorentz invariance was pointed out 
in [59] from the analysis of the zero temperature worldsheet Green's function. Notice that 
there is still a dependence of the widths on the modulus of the relative velocity between 
the heavy quarks and the plasma, implicit in T S f. 

In the non-conformal case, in addition to the Lorentz contraction, the velocity depen- 
dence of the drag coefficient affects to the distribution in the longitudinal direction, so 
it is not longer spherical in the rest frame of the quark distribution. This is even more 
explicit if the drag coefficient is taken as a function of the average longitudinal momentum 
Pl = "fvM. Then, the longitudinal momentum transfer has a rather simple form 

„3^ ..(, , P L 9fi 



"^'M' + fatJ' <448) 

A measure of the anisotropy in momentum space as a function of the velocity for different 
masses in the M = 2* theory can be found in figure |j. 



5. Jet quenching 

A different estimation of the energy loss of relativistic hadrons in the plasma [66] is based 
on a weak coupling calculation of the gluon radiation emitted by a dipole source in the 
limit E S> Pt S> T, where pt is the typical transverse momentum of the radiated gluons 
and E is the energy of the source. The momentum transfer to the medium q, known as 
'jet quenching' parameter, can be associated to a rectangular adjoint Wilson loop with two 
lightlike sides of length L_, related to partons moving at relativistic velocities, and two 
spatial sides of length L x <C L_ related to the transverse momentum of the radiated gluons 
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Figure 4: Anisotropy of the momentum distribution of heavy quarks in the rest 

frame of the ensemble as a function of the rapidity rj = arctanhv, for m/T = 
0, 1.28, 2.08, 3.42, 5.00, 6.94, 9.40, 12.58, 16.80 from bottom to top. 



Pt ~ 1/L X . In the planar limit, the value of the adjoint Wilson loop can be obtained from 
a Wilson loop in the fundamental representation 

(W Adj (C)) = \(W(C))\ 2 , (5.1) 

In [12] it was proposed to use the expectation value of the fundamental Wilson loop as a 
definition for the momentum transfer q at strong coupling 

(W(C)) = e ~^ LlL - , (5.2) 

and, using a classical string in the AdS^ black hole, the following result for the M = 4 
theory was found 

W=a= W s 4j ^T 3 . (5.3) 

One can see that once the dependence of the coupling and the temperature have been fac- 
tored out, the jet quenching differs from the momentum transfer obtained from worldsheet 



fluctuations of the trailing string ( 4.40 ). The relation between the jet quenching and the 
drag will be made explicit below. 



The expectation value of the Wilson loop (5.2) can be computed evaluating the classical 



action S(C) of a string ending on the contour C at the boundary. The contour C is chosen to 
be rectangular in the y~ = (t — y)/\^2 and x directions, with the length in the y~ direction 
very large, so it is possible to neglect any dependence on the y~ coordinate. The choice 
of coordinates for the associated string is r = y~ , r = a and the profile is determined by 
the function x(r). The Lagrangian is 

1„A 



C = ^^Vl - e 2B Vl + e 2A x' 2 . (5.4) 
a/2 
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From this one can derive the equation of motion 



ce 



^ 4 e«(l - e 2B ) - c 2 ' 
This can be simplified by defining a 'drag coefficient' function as 

rtr) = (l-e 2B ^) 1 / 2 n(r) 2 e 2A( - r K 



(5.5) 



(5.6) 



Notice that fi(r) at fixed r coincides with the drag coefficient ( f4.13| ) (up to the 2na' factor) 
at velocity v 2 = e 2B ^ r \ The radial coordinate is related to the energy scale in the dual 
theory, so there is a scale associated to each velocity. The map to a physical scale in the 
field theory could be ambiguous, depending on the choice of radial coordinate. A better 
option could be to use the energy scale defined by the average transverse kinetic energy 



of the heavy quark, using the momentum transfer ( 4.40 ) and assuming a non-relativistic 
dispersion relation as in 

(H) 



(E T ) 



Ml 
2M 



2^ 



St 



(5.7) 



The string profile ends at a finite value of the radial coordinate r > unless the 
integration constant c 2 < /i 2 (0) = /2q. This implies that in the non-conformal case there 
are two branches. In the c 2 > fi^ branch the action of the string grows linearly with the 
length L x in the x direction, while in the c 2 < /2q branch it grows quadratically. A similar 
behavior has been observed previously [67]. A numerical calculation in the J\f = 2* case 
shows that the branch with c 2 > /2q has a larger action for a given length, so it is disfavored. 
In the following only c 2 < /1q will be considered. 

Notice that x' ^ oo at the horizon. In principle this looks like the solution is irregular 
or should be continued beyond the horizon. However, it is possible to see that the solution 
corresponds to a regular string touching the horizon once Schwarzchild coordinates are 
used. The Gaussian radial coordinate r is related to the Schwarzschild radial coordinate u 
as dr = duj W g(u), where g(un) = at the horizon uh > 0. Then, 



dx I dx 
-dr = ^ g{u) Tu 



(5.8) 



If x' is constant at the horizon the derivative with respect to the Schwarzschild coordinate 

dx 
da 



diverges ^ — > oo, and this gives the right condition for a regular solution touching the 



horizon. 

With this choice, the length in the x direction and the regularized action of the string 

are 



T=y dr 



e~ A c 

\J[L 2 - C 2 



S 
2 



2^/2' 



Tta' Jo 



T —A ~ 

dr e fi 



(5- 
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Figure 5: Jet quenching in a s T units as a function of m/T. 



The jet quenching is obtained from the L x — > behavior, this corresponds to the c —* 
limit in the integrals. 



A 00 , 

— ~ c / ar 

2 7o At 



S 



c 2 L_ 



oo g -A 



2 4^/2^0' Jo 
So the (inverse of) the jet quenching coefficient is 



dr- 



ira 



dr- 



(5.10) 



(5.11) 



Taking into account the relation between the radial coordinate and the energy scale in 
the dual theory, the formula above could be interpreted as an average of the drag over 
all scales, showing a relation between the drag and the jet quenching. Using the relation 
between the scale and the velocity v = e B ^ , ( |5.11| ) could also be interpreted as an average 
over velocities, although the explicit form is more complicated. The value of q in units of 
the spatial string tension and temperature are plotted in figure |5| for different masses in 
the N = 2* theory. One can see that the jet quenching grows with m/T, the change is 
about ~ 50% for m/T ~ 17. The numerical value at m/T = q ~ 4.794 T a s is in good 
agreement with the analytical result (| 



20Fr (1) 
r(f) 



To* ~ 4.793 Ta, 



(5.12) 



6. Discussion of the results 

In [68] some experimental observables for heavy quarks were proposed such that it would 
be possible to make comparisons with perturbative QCD predictions. A matter of concern 
is how the comparison with 'real world' physics should be made. There are two parameters 
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that enter explicitly in both the drag and the jet quenching, the 't Hooft coupling and the 
temperature 7 . There have been different proposals to fix their value, for instance by fixing 
some physical quantity in both theories, like the energy density or the diffusion coefficient 
e.g. [8, 19]. However, different choices of the observables used to do the matching will lead 
to different quantitative results. Besides, it is difficult to quantify the value of m/T that 
would lead to the best comparison with QCD. A rough estimation of the deviation from 
conformal invariance that corresponds to QCD would be to compare the speed of sound 
or the equation of state of both theories 8 . The speed of sound is the most straightforward, 
although it would be desirable to use the equation of state as well to check its robustness. 
Using the values obtained in the lattice [3] and in the M = 2* theory [47], one finds that 
for T ~ 1.5 — 2 T c the value of m/T is in the approximate range 3 < m/T < 6. The 
uncertainty of this comparison motivates the study of quantities within the same theory, 
as has been done for the transport coefficients of the plasma, the shear to entropy density 
ratio being the most obvious example. 

The expression ( 4,16| ) shows that the ultrarelativistic drag coefficient is sensitive to 



the infrared physics, and that it should be related to transport coefficients of the plasma 
obtained from the energy momentum tensor. One can also see from the formulas (|3.6|) and 



( 4.18 ) that a smaller spatial string tension in the non-conformal theory implies a relatively 
larger drag at high velocities and vice versa. A relation between the jet quenching and the 
drag coefficient is given in section [f| in the formula ( 5.11 ) one can see that the value of the 



jet quenching is an average of the drag coefficient over all scales. Therefore, their values 
are correlated as should be expected. 

The numerical results for the Af = 2* theory presented in figures |2| and [| show that 
the energy loss is larger in string tension units than in the non-conformal case, and grows 
as the mass deformation becomes more important. For the range of masses studied here 
< m/T < 17, the jet quenching changes up to ~ 50%. The drag coefficient can be more 
than twice as large as the conformal value at high velocities. 

In addition to the drag coefficient, the momentum broadening of heavy quarks [59,60] 
was studied in section The interpretation of the original results is much less clear in this 
case, since the values that were obtained did not fit with the expectations from a Langevin 
description of the dynamics of heavy quarks at thermal equilibrium with the plasma. A 
recent work [63] , gives a possible interpretation in terms of energy loss through non-thermal 
radiation processes. The transverse and longitudinal momentum transfer coefficients are 
found to be 

K T = 2 7 T stM , KL = 2 7 3 T^l + ^^0 . (6.1) 

where Pl is the average longitudinal momentum of the heavy quark. It is argued that these 
formulas are compatible with steady state configurations of quarks subject to an external 
force and moving through a relativistic viscous medium. The additional 7 2 factor in kl 
relative to kt is argued to be due to Lorentz contraction in the longitudinal direction, 



7 In the Af = 2* theory the mass deformation is an additional parameter that does not enter explicitly 
but that modifies the value. 

8 I would like to thank Krishna Rajagopal and Hong Liu for suggesting these comparisons. 
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under the assumption that the dynamics of heavy quarks can be described with a Liouville 
equation, at least as a first approximation. 

In the conformal case this implies that the distribution in momentum will be spherical 
in the rest frame of the quark ensemble, but in the non-conformal case there is an anisotropy. 
It is surprising that the anisotropy of the system does not introduce an anisotropy in the 
distribution in the conformal case. It would be interesting to understand whether this 
can be derived from conformal invariance or if it is a characteristic of the holographic 
computation. It is also possible that, even in a first approximation, the dynamics of 
the heavy quarks are described by something more complicated than a simple Liouville 
equation, in such a way that the simple relation of the 7 2 factor to a Lorentz transformation 
does not hold. 

In the J\f = 2* theory a broadening of the longitudinal momentum is observed (Fig. Q), 
related to the growth of the drag coefficient at larger velocities. A broadening of the 
momentum distribution in the longitudinal direction has also been observed in perturbative 
calculations of quarks moving through an anisotropic plasma [69]. This seems to be the 
right qualitative behavior that is consistent with the observed distribution of jets in the 
plasma. It would be interesting to do a comparison with the broadening predicted by 
holographic models. 
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A. Numerical method 



The numerical calculations are made following the conventions of [46,47]. The Einstein 
equation ( 2.1l| ) together with the equations of motion for the scalar fields ( [2. S| ) form a 
system of second order ODEs with four independent equations. The radius of curvature is 
fixed to L = 2 to make the gauged supergravity coupling unity. The leading terms in the 
asymptotic expansion of the scalar fields as r — > oo are 



— r 11 



a 

73 



—f= — -Poore + 



From (|2,13| ) , solutions with Af = 2 supersymmetry should satisfy the condition 



Xoo 



3p 







(A.l) 



(A.2) 



The asymptotic behavior of the metric functions should be the same for all solutions. From 
the holographic dual perspective this means that they flow to the same ultraviolet fixed 
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point. Comparing to the usual AdS^ black hole, as r — ► oo, 

e A(r)^ e r/2^ g S(r) _^ ^ (A _ 3) 

This fixes completely the boundary conditions at infinity for a fixed value of the mass. 

Regularity of the solution at the horizon imposes additional constraints as r — > 0, 
namely 

e B(r) ^ 5 Qr + . . . ; A /( ) = o, «'(0) = 0, X '(0) = . (A.4) 

It will be convenient to extract the non-analytic behavior of the function B(r) as 

B(r) =4b(r) - AA{r) + log(6 r), 6'(0)=0. (A.5) 

It is actually easier to implement these conditions at the horizon and shoot toward the 
boundary, this is the approach used here. The conditions ( A.3 ) can be imposed by shifting 
the functions A(r) and B(r) afterwards. The condition (|A.2|) is used as a criterion to keep 
the solution as the space of initial values x(0) = Xo and a(0) = «o is explored 9 . In the 
numerical solutions, the value of the supersymmetric condition is always kept under 10 -2 . 
With this condition, the numerical solutions have a small relative mass difference between 
fermions and bosons 

' 771,5 ~ m f 



< 0.005. (A.6) 

m b 

The shooting is done in Mathematica using a set of six first order differential equations 
and a second order differential equation. The choice is inspired by the equations of the su- 
persymmetric flow (^1]) . Following those, the deviation from extremality is parameterized 
by three functions a(r), c(r) and b\(r) that vanish as the non-extremal flow approaches the 
supersymmetric one when r — > oo, 

, IdW , IdW w 

a =4^ + a ' X = 4fe +C ' 6 = Y + &1 ' (A - 7) 

where W is the superpotential (p.7|). The Einstein equations and the equations of motion 
for the scalar fields give 

/ ( -\ At. 4 rxr 1 d 2 W\ IdW , dW 1 d 2 W n 

ra +a { 1+ Arh - r w + J + 4^ + rh ^ + rc 4 = ' (A - 8) 

/ (, , , 4 1 d 2 W\ IdW , dW 1 d 2 W n , An , 

rc +c { 1+ 4rh - r w + 4 r az ) + -^ +rb ^ +ra 4o^=^ (A - 9) 

/ 8 \ 2 1 dW 1 dW 

rb[ + h 2 + 4r6i - -rW - -W - -rc- -ra— = . (A.10) 

\ 3/3 3 ox 3 da 

Finally, there is a second order equation that can be solved once the solutions to the first 
order system are known 

rA" + A'- 2b' + ArA'b' + 4rb' 2 - rb\ +-r( + ^-o/] = . (A.1I) 

3 \dx da J 



9 1 would like to thank Steve Paik for providing a set of initial values that helped very much for this 
analysis. 
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The boundary values at r = of the extra functions are fixed by the requirement of 
regularity 



a(0) = - 



1 dW 



4 da 



, c(0) = - 



1 dW 



a=a ,x=XO 



4 d x 



6i(0) = 



a=a ,X=X0 



. (A.12) 



a=a ,X=XO 



There are only two more boundary conditions that need to be fixed 

6(0) =0, A(0) = 0. (A.13) 
This choice implies that the solutions will correspond to different temperatures. 
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